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FOREWORD 


This  Special  Report  represents  the  work  performed  by  Thiokol  Chemical 
Corporation,  Wasatch  Division,  Brigham  City,  Utah,  under  Contract  F04611-72- 
C-0015  from  1  Nov  1971  thru  29  Feb  1972. 

Technical  direction  was  provided  by  the  Air  Force  Rocket  Propulsion  * 

Laboratory  at  Edwards  Air  Force  Base,  Edwards,  California.  * 

t 

This  technical  report  has  been  reviewed  and  is  approved.  * 

R.  G.  Schuder,  Capt/USAF 
Project  Engineer/AFRPL 


ABSTRACT 


This  report  presents  a  theoretical  stress  analysis  of  OOAMA/RPL 
motor  S/N  0012257.  The  analysis  involved  motor  tip-over  from  horizontal 
to  vertical  position  and  rerun  of  the  board  course  reported  previously  in 
Thiokol  document  TWR-4728.  The  results  reported  represent  a  major  step 
in  the  use  of  instrumented  motors  to  verify  the  accuracy  of  modulus  depend 
ent  predictions  of  rocket  motor  grain  structural  behavior. 
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I.  INTRODUCTION  AND  SUMMARY 


This  report  presents  the  results  of  a  theoretical  stress  analysis  of  the 
OOAMA/RPL  stress  motor,  S/N  0012257.  The  data  presented  forma  continuation 
of  the  analysis  described  in  ref  1. 

The  test  conditions  considered  in  this  analysis  were: 

1.  Tip-over  from  horizontal  to  vertical  position 
and  back  to  the  original  horizontal  position.  The 
motor  remained  in  the  vertical  position  for 
approximately  3  days. 

2.  Rerun  of  the  board  course  test  as  described  in 
ref  1. 

The  tip-over  analysis  includes  the  effects  of  the  change  from  horizontal  to 
vertical  slump  loading.  In  addition,  whenever  possible,  an  attempt  was  made  to 
compensate  for  thermal  gradients  that  occurred  within  the  motor  during  the  tip- 
over  test. 

The  predicted  variations  in  the  stresses  and  strains  experienced  by  the 
motor  during  the  tip-over  from  horizontal  to  vertical  are  dependent  on  grain 
location.  The  maximum  predicted  variations  were  as  follows: 


1. 

Normal  stress  at  the  case  (psi) 

2.59 

2. 

Shear  stress  at  the  case  (psi) 

0.87 

3. 

Bore  strain  {%) 

2.71 

In  addition,  the  following  variations  were  predicted  due  to  thermal  gradient 

changes: 

1. 

Normal  stress  (psi) 

-0.4 

2. 

Bore  strain  (%) 

-1.23 

No  attempt  was  made  to  predict  thermal  gradient  effects  on  shear  stress.  Table  I 
shows  the  results  of  the  tip-over  analysis  for  various  times  during  the  test  for  the 
applicable  gages. 

The  dynamic  analysis  was  based  on  the  same  analysis  techniques  used  in 
the  previous  analysis  (ref  1).  Various  loading  conditions  were  obtained  by  travers¬ 
ing  a  spaced  board  course  with  a  GFE  transporter  pulled  by  a  tractor.  Vertical 
accelerations  (single  amplitude)  measured  in  the  propellant  grain  varied  from 
0.31  to  1. 11  g*s.  Resultant  predicted  values  of  stress  and  strain  are  shown  in 
Table  II.  Additional  accelerometers  mounted  on  the  starpoints  indicated  that 
the  entire  grain  responds  to  the  same  transportation  frequency.  No  measurable 
phase  lag  was  observed  in  any  of  the  accelerometers. 


II.  ANALYSIS  METHOD 

Budget  and  time  considerations  dictated  that  large  computer  runs  be  kept 
to  a  minimum.  Previously  generated  grids  and  runs  were  used  wherever  possible. 
However,  an  axisymmetric  equivalent  port  area  grid  (Figure  1)  was  required. 
Three  parameters  are  of  interest  in  this  analysis. 

1.  Normal  stress  at  the  bondline  (crr),  gages  Nl, 

N2,  and  N3  located  at  planes  70,  31,  and  70, 
respectively. 

2.  Shear  stress  (T),  gages  SCI,  SC2,  SI,  S2,  and 
S3  located  at  planes  70,  31,  200,  60,  and  200, 
respectively. 

3.  Bore  strain  at  the  star  tip  (t  g  )>  gages  Cl  thru 
C8,  located  at  planes  31  and  70. 

Gage  .locations  are  shown  in  Figures  2  and  3. 
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Computer  programs  utilized  in  the  analysis  were: 

1.  S3016,  A  Finite  Element  Plane  Strain  Stress 
Analysis  Program,  ref  2. 

2.  S3036,  An  Axisymmetric  Finite  Element  Stress 
Analysis  Program,  ref  3. 

3.  S3053,  Elastic  Hollow  Cylinder  Stress  Analysis 
Program. 

4.  S3069,  Linear  Viscoelastic  Hollow  Cylinder  Pro¬ 
gram,  based  on  theory  presented  in  ref  4. 

A.  TIP-OVER  TEST  ANALYSIS 

1.  Horizontal  Slump  Load 

The  analysis  assumes  that  all  stresses  and  strains  are  zero  when  the  motor 
is  in  the  horizontal  position  (gages  calibrated  to  read  zero).  Actually  the  bulk  of 
the  grain  is  responding  to  a  1  g  horizontal  slump  load.  When  the  motor  is  positioned 
vertically,  thelg  horizontal  slump  load  will  relax  or  subtract  from  the  zero  gage 
readings.  This  effect  must  be  taken  into  account  in  analysis  of  the  vertical  slump 
load. 

This  analysis  was  made  with  the  use  of  Program  S3016.  A  1  g  slump  load 
was  applied  to  the  propellant  grain  with  a  body  force  applied  uniformly  around  the 
case  to  obtain  an  equilibrium  force  balance.  A  180  deg  symmetrical  sector  of  a 

t 

cross-sectional  slice  through  the  motor  was  utilized.  The  resulting  input  grid  is 

shown  in  Figure  4.  Parameters  obtained  from  this  analysis  were  a  ,  ,  €  . ,  and  T 
)  .  rh  0h 

predicted  -fop  each  applicable  gage  location. 

2.  Vertical  Slump  Load 

This  analysis  consisted  of  a  run  on  Program  S3036  utilizing  an  axisymmetric 
grid  with  an  equivalent  port  area  smooth  bore  geometry  (Figure  1).  Loading 
consisted  of  a  1  g  vertical  slump  reacted  at  the  aft  skirt  and  a  1,500  lb  load 
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applied  at  the  forward  skirt  to  incorporate  the  effects  of  the  forward  lifting  ring. 

Parameters  obtained  from  this  run  were  a  ,  €  ,  and  T  for  each  applicable 

rv  va  v 

instrumentation  plane.  Two  additional  runs  were  made  with  a  thermal  load  applied. 

These  were  a  plane  strain  run  with  Program  S3016,  with  the  geometry  and  grid 

shown  in  Figure  5  and  a  plane  strain  equivalent  port  area  cylinder.  Parameters 

obtained  from  these  runs  are  c  and  €  respectively. 

0  s  0  c 


These  data  were  then  used  to  obtain  a  concentration  factor  (K)  for  subse¬ 
quent  adjustments  of  the  axisymmetric  output  to  the  star  geometry.  These  adjust¬ 
ments  were  made  for  each  plane  of  interest  as  follows: 


c 


ev 


K  € 


0va 


where  K  =  e  /( 

0s  0c 


3.  Thermal  Gradient  Analysis 

The  motor  experienced  significant  thermal  gradient  changes  during  the 
tip-over  test.  Outside  temperatures  dropped  to  a  low  level  during  the  test  causing 
the  case  and  propellant  near  the  case  to  drop  to  a  temperature  considerably  lower 
than  the  storage  temperature. 


Normal  stress  and  bore  strain  data  were  adjusted  to  reflect  these  temperature 
changes  at  several  times  during  the  test  cycle. 

The  following  runs  were  required  to  make  the  above  described  adjustments. 


1.  Axisymmetric  equivalent  port  area  (Figure  1) 

with  uniform  thermal  load.  The  output  parameters 

were  a  and  € 

rax  0ax 

2.  Plane  strain  star  geometry  (Figure  6)  with  uni¬ 
form  thermal  load.  The  output  parameters  were 

o  and  € 
rs  0s 

S.  Plane  strain  equivalent  port  area  cylinder  with 

uniform  thermal  load  identical  to  a  and  b.  The 

output  parameters  were  <7  and  €  • 

rc  0c 
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4.  Plane  strain  cylinder  geometry  with  equivalent  port 
area.  These  runs  used  temperature  gradient  data 
for  each  time  considered  and  were  obtained  from 
thermocouples  T9  thru  T13  (Figure  3).  The  output 
parameters  were  c t  ,  C  . 

rt  0t 

From  these  runs  two  concentration  factors  are  calculated  for  each  parameter 
to  be  adjusted. 


T.,*,se  concentration  factors  are  then  used  to  obtain  the  adjusted  stress  and 
strain  values  due  to  thermal  gradients 


=  K 

K  , 

a 

rtg 

ag 

at 

rax 

0tg 

=  K 

eg 

Ktt 

C  6ax 

The  values  obtained  from  the  horizontal  slump,  vertical  slump,  and  thermal 
gradient  analyses  were  then  combined  to  obtain  the  predicted  delta  values  of  cr^>  e  g, 
and  T  at  the  selected  times  and  gage  locations.  All  stress/ strain  values  were 
considered  0.0  prior  to  lift  to  the  vertical  position. 


Material  Properties 


Material  properties  used  for  this  analysis  (see  Table  HI)  were  the  same  as  those 
used  in  the  previous  analysis  (ref  1)  except  for  the  equivalent  viscoelastic  propellant 
modulus.  Also  the  Thiokol  value  for  thermal  coefficient  of  linear  expansion  (2) 
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was  used.  This  value  was  chosen  because  it  resulted  in  better  correlation  with  test 
data  in  the  previous  analysis. 

A  procedure  has  been  developed  to  reduce  the  required  viscoelastic  analysis 
to  an  equivalent  elastic  analysis.  This  is  accomplished  by  the  use  of  Programs 
S3069  and  S3053.  Input  geometry  to  both  programs  consists  of  a  cylinder  with  a 
15.2  in.  inner  radius,  a  32.7  in.  outer  radius,  and  a  0.144  in.  case  thickness. 

A  viscoelastic  analysis  was  run  on  Program  S3069.  Input  to  this  program  consists 
of  relaxation  modulus  (E  ^)  versus  temperature  reduced  time  (see  Figure  6)  tempera¬ 
ture  shift  factor  (A^)  versus  temperature  (see  Figure  7)  and  time-temperature  history 
of  the  motor.  Output  from  this  program  is  cx^  at  the  case/propellant  interface  as 
a  function  of  time.  An  elastic  analysis  was  then  run  on  Program  S3053  for  thermal 
loads  experienced  during  the  test  and  various  values  of  modulus.  Output  from 

these  runs  gives  cr  as  a  function  of  modulus.  The  output  values  from  these  two 
r 

programs  are  then  combined  to  obtain  a  modulus  (equivalent  elastic  modulus) 
which  will  give  the  same  interface  stress  in  the  elastic  cylinder  as  is  obtained  from 
the  viscoelastic  cylinder  analysis.  Due  to  the  complex  thermal  loading  during  the 
test,  an  effort  was  made  to  obtain  an  average  value  of  E  during  the  test  cycle.  This 
value  is  160  psi  and  was  used  for  all  calculations  during  the  test. 

B.  DYNAMIC  ANALYSIS 

Vibration  loads  were  induced  in  the  OOAMA/RPL  stress  motor  by  travers¬ 
ing  a  spaced  board  course  with  a  GFE  transporter  pulled  by  a  tractor.  The  load 
was  varied  by  making  a  number  of  runs  at  different  speeds. 

The  motor  case  had  three  triaxial  accelerometers  (see  Figure  3):  one  at  each 
end  of  the  motor  and  the  third  located  50  in.  aft  of  the  head  end  (plane  50).  A 
vertical  accelerometer  was  mounted  on  the  propellant  at  plane  50  (V4)  15  deg  off 
the  motor  vertical  axis.  This  accelerometer  was  used  to  determine  the  loads  used 
for  the  analysis. 
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The  previous  analysis  (ref  1)  showed  that  predicted  stress  and  strain  values 
are  linear  functions  of  modulus  (E)  and  acceleration  loads  (g).  The  relationships 
derived  from  these  data  are: 

<Tr2  =  arl  (62/gl) 

T2  =  T1  (g2/gl) 

=  eei  (g2/gl)  (E1/E2) 

Where  the  subscript  1  refers  to  the  reference  parameters  and  the  subscript  2 
refers  to  the  new  parameters. 

These  equations  in  conjunction  with  the  modulus  shown  in  Table  II  were 
used  to  obtain  the  results  also  shown  in  Table  II.  The  data  obtained  from  ref  1  was: 

a  =  1.56  psi 

r 

T  =0.66  psi 

€  =  0.07  percent 

0 

g  =  0.8  g 

E  =1, 860  psi 

Dynamic  modulus  (E1)  was  obtained  from  Figure  4  where  E*  equals  3G*.  The  curve 
of  G*  versus  frequency  (Figure  8)  based  on  a  JANNAF  modulus  of  800  psi  was 
used  whereas  the  previous  analysis  used  the  G'  curve  based  on  a  JANNAF  modulus 
of  500  psi  (page  19,  ref  1).  This  should  give  a  more  realistic  value.  Other 
material  properties  used  in  this  analysis  were  the  same  as  the  previous  analysis 
and  are  shown  in  Table  HI. 


ID.  ANALYSIS  RESULTS 


A.  TIP-OVER  ANALYSIS 

The  results  of  this  analysis  are  summarized  in  Table  I.  Values  are  reported 
for  the  applicable  gages  and  for  different  loading  conditions  which  occurred  during 
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the  test.  The  stresses  and  strains  shown  are  predicted  delta  values  from  the  time 
just  prior  to  lifting  the  motor  to  the  vertical  position  (5  Nov  1971  at  1745  hours). 

The  final  values  shown  are  for  the  time  (8  Nov  1971  at  1415  hours)  just  after  the 
motor  was  lowered  to  the  initial  horizontal  position.  It  should  be  noted  that  the 
final  values  did  not  return  to  a  0  delta  value  due  to  thermal  changes  that  occurred 
in  the  motor  during  the  test. 

B.  DYNAMIC  ANALYSIS 

This  analysis  was  performed  for  selected  times  during  each  run  over  the 
board  course  (two  times  from  Runs  1  and  5  and  1  time  for  Runs  2  and  3).  These 
results  are  summarized  in  Table  II  and  reflect  only  the  stress  and  strain  components 
due  to  vertical  acceleration  loading.  Measured  case  longitudinal  and  tangential 
strains  are  shown  in  Table  IV. 


IV.  CONCLUSIONS 

The  analytical  prediction  of  stresses  and  strains  to  be  measured  in  the 
OOAMA/RPL  motor  program  have  been  presented.  Although  the  analysis  was  not 
very  refined,  the  predictions  represent  a  reasonable  best  estimate  of  the  stresses 
and  strains  to  be  measured. 

The  significance  of  all  the  OOAMA/RPL  stress  motor  programs  is  that 
they  will  provide  grain  stress  analysts  with  information  previously  unattainable. 

In  the  past,  the  grain  analysts  could  not  verify  the  accuracy  of  modulus  dependent 
predictions  because  no  means  of  verifying  these  predictions  existed.  The  advent 
of  usable  stress/strain  gages,  however,  now  makes  such  means  possible.  Future 
investigations  of  instrumented  motors  will  significantly  extend  the  grain  analysts' 
knowledge  of  grain  structural  behavior.  The  program  presented  in  this  report 
represents  a  major  step  towards  this  achievement. 
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Figure  2.  Propellant  Stress,  Strain  and  Temperature  Instrumentation  Location 


Figure  5.  Thermal  Stress  Analysis  Grid 


LOG  TIME,  LOG  SHIFT  FACTOR,  LOG  t/AT  (MIN) 


Figure  8.  Dynamic  Shear  Modulus  vs  Frequency 
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TIP-OVER  TEST  ANALYTICAL  STRESS  ANALYSIS  RESULTS 


BOARD  COURSE  VERTICAL  RESPONSE  AND  STRESS  ANALYSIS  RESULTS 


TABLE  III 


Tip-Over  Analysis 
Poisson's  Ratio 
Modulus  (psi) 

TCLE  (a)  (in./in. /°F) 
Dynamic  Analysis 
Poisson's  Ratio 


MATERIAL  PROPERTIES 


Case 

Liner 

Propellant 

0.3 

0.499 

0.499 

6 

30.0  x  10 

130 

160 

p 

7.076  x  10 

1.5  x 10“4 

54.0  x  10H 

0.3 

0.499 

0.499 

6 

30.0  x  10 


Modulus 


480 


(See  Table  II) 


BOARD  COURSE  LONGITUDINAL  AND  TANGENTIAL  CASE  STRAIN 
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